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(57) Abstract: A method of creating an image of brain fibers 
includes exposing tlie brain fibers to a magnetic resonance imag- 
ing process. The data acquisition from the magnetic resonance 
imaging includes the acquisition of diffusion-weighted images 
that are later employed to calculate an apparent diffusion con- 
stant at each pixel along more than six axes. The data is in- 
troduced into a microprocessor which calculates six variable in 
a diffusion tensor and obtains a plurality of eigen values and 
eigen vectors. This may be accomplished by employing a dif- 
fusion sensor which is diagonalized to obtain three eigen values 
and three eigen vectors with the six values being subjected to 
further microprocessing to generate imaging information repre- 
senting the properties of the fibers. The process in a preferred 
embodiment includes the initiation of fiber tracking by select- 
ing a pixel for initiation of the same, connecting of pixels and 
effecting a judgement regarding termination of the pixel track- 
ing in each direction based upon the randomness of the fiber 
orientation of the adjacent pixels. 
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METHOD OF FIBER RECONSTRUCTION EMPLOYING 
DATA ACOUIRED BY MAGNETIC RESONANCE IMAGING 

BACKGROUND OF THE INVENTION 

5 

1. Field of the Invention 

The present invention provides a method of employing data obtained 
through magnetic resonance imaging in creating three-dimensional brain fiber 
reconstructions and, more specifically, it provides such a method which is particularly 
10 useful in respect of white matter fibers. 

2. Description of the Prior Art 

It has been known for many purposes to attempt to image portions of a 
human or animal's brain for diagnostic, research, or there^utic purposes. In order to 
understand the activity of and relationships between multiple cortical regions of the 

15 brain, it is important to analyze their physical connectivities, such as axonal fibers. 
The projections of these fibers can be traced in experimental animals by observing 
axonal degeneration following carefully placed experimental brain lesions or, 
alternatively, by injecting and subsequentiy localizing radio-isotopes or other chemicals 
which are taken up by nerve cells and actively transported along their axons. 

20 Comparable human data are necessarily much more limited because they have been 
obtained only after postmortem examinations of patients with naturally occurring 
lesions such as iiyuries or infarcts. 

As a result of the foregoing, progress in understanding the structure of 
cognitive association pathways and their dysfunction in many disorders has been slow 

25 due to the lack of a non-invasive method for fiber tracking. Conventional MR imaging 
can distuiguish gray and white matter, but provides no information about the fiber 
orientation in white matter. Diffusion-weighted MRI allows in vivo mapping of the 
anisotropic and isotropic diffusional properties of brain water, and has revealed a high 
degree of diffusional anisotropy in white matter. Although this finding has been 
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teixtatively attributed to preferential water diffusion along axons and/or flieir myelin 
sheaths, it has not been known to show tracking of neuronal projections. 

There remains, therefore, a need for improved methods of imaging 
white matter fibers in the brain including the brain of a living human being. 

SUMMARY OF THR INVENTION 

The present invention involves a method of creating an image of brain 
fibers which includes the data acquisition by the magnetic resonance imaging process 
and the data processing to generate imaging ioformation relating to the fibers. 

The data acquisition by the magnetic resoi}ance imaging contains the 
acquisition of so-called diffusion-weighted images that are later used for the calculation 
of apparent diffusion constant at each picture element (pixel) along more than six axes. 
This can be accomplished by using a pair of magnetic field gradients to sensitize the 
magnetic resonance imaging. Conventional magnetic resonance imaging scaxmers are 
equipped with three magnetic field gradient units. By combining these units and by 
changing the strength, a series of images sensitized to water diffusion along desired 
direction can be recorded. 

After the data acquisition, the imaging information is transferred to the 
conq)uter for the fiber analysis. First, intensity of eadi pbcel of diffiision-weigjited 
images with various gradient combination and strength are fitted to calculate six 
independent variables in a 3 x 3 diffusion tensor. The dif^ion tensor is then 
diagonalized to obtain three eigen values and three eigen vectors. These six values are 
subjected to the further computer processing to generate images representing the 
properties of the fibers. This processing preferably consists of three parts: initiation of 
brain fiber tracking; pixel connecting; and the judg^nent of the termination of the 
fibers. For example, a tracking of projections of fibers can be initiated from a point in 
a three-dimensional space arbitrarily chosen by a user and propagated in both 
directions according to the direction of the fiber (the eigen vector associated with the 
largest eigen values). Each time the tracking leaves a pixel to the next pixel, 
judgement is made whether the fiber is continuous or terminated based on randomness 
of the fiber orientation of the adjacent pixels. 
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It is an object of the present invention to provide a method for using 
data acquired by magnetic resonance imaging of a brain in creating a three-dimensional 
fiber structure for white matter fibers of the brain. 

It is a further object of the present invention to provide such a method 
wherein axonal fiber may be imaged so as to provide images containing details 
regarding the patient white matter fibers. 

It is a further object of the present invention to provide an automated 
means of converting magnetic resonance imagmg data into bram fib^ three- 
dimensional image. 

It is a further object of the present invention to provide such a method 
which may be practiced on living human beings. 

These and other objects of the invention will be more fully understood 
from the following detailed descriptions of die preferred embodiments on reference to 
the illustrations appended hereto. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figures la and lb illustrate principles of water diffusion anisotropy as 
known in the prior art. 

Figure 2 illustrates a flow diagram showing a sequence of operations of 
the present invention. 

Figures 3a through 3c are schematic diagrams of fiber tracking by a 
first method of the present invention. 

Figures 4a and 4b illustrate, respectively, two and three-dimensional 
presentations of the parietal lobe of a rat brain with representations of prominent brain 
fiber bimdles. 

Figures Sa through 3f illustrate three-dimensional projections and two- 
dimensional validation of 8 fiber bundles in a rat brain. 

Figure 6 is a schematic diagram of fiber tracking by a second method of 
the present invention. 

Figure 7 is an example of an optic tract reconstructed by the second 
method of this invention. 
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DESCRIPnON OF THE PREFERRED EMBODIMENT 
As used herein, the term ""patient" refers to members of the animal 
kingdom including human beings. 

The terms "fiber" or "fibers" as used herein shall refer to a group of 
5 axonal fibers which have at least portions thereof projecting in generally the same 
direction. 

The principle of water diffusion anisotropy of the prior art as compared 
with axonal projections of the present invention is shown in Figure 1 . Figure la is a 
schematic view of restricted water diffusion (soUd sphere) in an environment with 

10 strongly aligned fibers (depicted by bars). The diffiision properties can be quantified 
by an ellipsoid (Figure lb) with three principal axes (Xi, X^, and X3) and a certain 
orientation of the main axis desoribed by a vector. For a region where axons are 
aligned (depicted by long cylinders in Figure la), water would be restricted in the 
direction perpendicular to the axons and consequentially diffuse preferentially in a 

15 direction parallel to them. The properties of such water diffusion can be represented : 
mathematically by an ellipsoid, as illustrated in Figure lb. This so-called diffusion 
ellipsoid can be fully characterized by the diffusion constants X2, and along the 
three orthogonal principal directions corresponding to the longest, intermediate, and 
shortest axes, respectively, and the three-dimensional (vector) direction of the longest 

20 axis. For example, a case for which > > X2 = A3 (anisotropic diffusion) suggests 
the existence of dense and aligned cylindrical axonal structures, whereas a result of Xx 
= ^2 = ^3 (isotrc^ic diffusion) suggests sparse or unaligned axons. When diffusion is 
anisotropic, the direction of indicates the preferential fiber orientation within that 
voxel. All ellipsoid parameters can be obtained from diffusion MRI measurements for 

25 each volimie element (voxel) of a bram image. The diffusion ellipsoid for a voxel 

represents an average over all neuronal structures contained within. Although 

anisotropics have been detected previously in white matter, it previously has not been 

known to relate these to neuronal trajectories. 

In order to accomplish the three-dimensional reconstmction of axonal 

30 fibers, the present invention employs a computer which tracks the fibers in such a 

4 
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three^limensioiial space. Figure 2 illustrates a flow diagram for the reconstruction 
process. The first block indicates data acquisition by magnetic resonance imaging 
which involves the acquisition of so-called diffusion-weighted unages that are later 
employed to calculate the apparent diffusion constant at each picture element (pixel) 
5 along more than six axes. This can be accomplished by usmg a pair of magnetic field 
gradients to sensitize the magnetic resonance imaging. Conventional magnetic 
resonance imaging scanners are equipped with three magnetic field gradient units. By 
combining these imits and by altering the strength, a series of images sensitized to 
water diffusion along the desired direction can be recorded. The image information is 

10 then introduced into a computer for the fiber analysis. In the diffusion tensor 
calculation form, the intensity of each pixel of diffusion-weighted images v^th various 
gradient conlbinations and strengths are fitted to calculate six independmt variables in 
a 3 X 3 diffusion tensor. The diffusion tensor is then diagonalized to obtain three eigen 
values and three eigen vectors. The images representing the properties of the fibers 

15 are then determined by further computer processing of the six values. In the fiber 
tracking portion of the method, as shown schematically in Figure 2, the input of an 
initiation pixel begins the tracking process which involves pixel connecting and 
ultimately judgement or a decision based upon the tracking reaching termination of the 
fiber in each direction. If the comparison in the computer results in the conclusion that 

20 the end of the fiber has been reached ("yes"), that is the end of the process as to that 
axonal fiber. If not ("no"), the fiber tracking process contioues. 

By way of an example, a tracking of projections of fibers can be 
initiated from a point in three-dimensional space which is arbitrarily selected by a user 
and propagated in both directions according to the direction of the fiber (the eigen 

25 vector associated with the largest eigen vahie). Each time the tracking leaves a pixel 
and proceeds to the next pixel, a judgement is made as to whether the fiber is 
continuous or termimted based on randomness of the fiber orientation of the adjacent 
pixels. 

In a first method of fiber reconstruction of the invention, with reference 
30 to Figure 3, the initial point is input by a user in an arbitrary manner and then tracking 

5 
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is performed in both forward and backward directions. The actual tracking consists of 
two stages. The first stage is a decision concerning the pixel connection and then the 
start and end points of projections are determined as by the flow chart of Figure 2. 
Figure 3 a shows an example of the method for the decision-making with respect to the 
5 pixel connection in which the tracking follows the direction of eigen vectors associated 
with the largest eigen values. The curve lines in Figure 3a represent a method of pixel 
coimection in whidi curved lines represent axonal projections and the small arrows the 
direction of eigen vectors associated with the largest eigen values of each pixel. 
Starting from a pixel indicated by an asterisk in Figure 3a, connected pixels are shown 

10 by shaded pixels. The endpoint of the projection is determined as shown in Figures 3b 
and 3c on the basis of the occurrence of a sudden transition of fiber orientation. The 
tracking is terminated when the environment is random. 

While the endpoint may be determined by several means, one approach 
to quantifying the severity of such a transition is to use the sxmimation of inner 

15 products of nearby data points. A three-dimensional axonal projection is tracked by ^ 
referencing the distance-weighted group orientation of nearby vectors Figure 3b. 
When vector orientation becomes random, as judged quantitatively from summation of 
the iimer products of these vectors, as in Equation 1, the tracking is ended Figure 3b. 
The computerized means of performing the process accomplishes this in two stages. 

20 In the present method, first, a decision concerning the axonal orientation is made. 
Then the start and end points of projections are determined. As shown in Figure 3b, 
MRI information on fiber orientation is discreet and each data point represents the 
average orientation of axons within a voxel. The fiber tracking, on the other hand, is 
performed in a continuous number field by referencing nearby three-dimensional 

25 discreet data grids. The endpoint of the projection is judged based on the occurrence 

* 

of sudden transitions of fiber orientations (Figure 3c). The severity of such a transition 
is quantified through a parameter R, presenting the summation of the inner products of 
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nearby data points: 

where v^i is the unit vector representing the longest principal diffusion axis Xj, the 
index i is the three-dimensional coordinate of a pixel of interest (fiber coordinate) and j 
5 are the coordinate of pixels surrounding the pixel i, s is the number of data points 
referenced. In general, the fiber coordinate, i, is considered to have 26 surrounding 
points, j. A point such as I or j may be identified by three numbers each related to 
reference axis. The term "abs" refers to the number having an absolute value without 
reference to " + " or 

10 In the experimental work, R values were calculated from the four 

closest data points and a fiber was judged to be discontinued for R values less than 0.8. 
There are ways to modify this criterion. For example, the number of data points 
referred and the method to weight the each point may be modified. The optimal way 
also depends on the relative size of the image resolution and brain structures. Because 

15 of this high threshold used (0.8), only fibers in a strongly aligned environment were 
tracked. The procedure of mappiog the neuronal connections is started through the 
input of an arbitrary point in three-dimensional space, i.e., "i," after which the extent 
and distribution of the axonal projections into functional regions is traced in both the 
orthograde (forward) and r^ograde (backward) directions as described above. As an 

20 example of this method, eight different pathways were traced in a formalm-fixed rat 
brain, as illustrated in Figures 4 and 5. With reference to Figure 4, two-dimensional 
and three-dimensional tracking of prominent axonal projections in the context of the 
present invention will be considered. As shown in Figure 4a, two-dimensional vector 
field presentation in the parietal lobe of a rat brain as localized firom a section of the 

25 Tj-weighted MR image. Regions for prominent fiber bundles are indicated by color 
lines (yellow : corpus caUosum (cc) and external capsule (ec), green : fimbria (fi), and 
red : internal capsule (ic). The presence of preferential axonal orientation in the gray 

7 
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matter which can be appreciated in Figure 4a suggests that Ihe fiber reconstruction 
scheme can also be applied to gray matter. 

Alternatively, the endpoint may is defined when the extent of 
anisotropy is weaker thaa a threshold value. The extent of the anisotropy can be 
5 quantified in various ways using three eigen values, A^, Xj, and A^, such as by 
calculadng the ratio between Xi and X^. 

In Figure 4b, there is shown a three-dimensional presentation of the 
fibers. Color-coding is the same as in Figure 4a except for the blue color which shows 
axons tracked fi*om the corpus callosum into the external capsule. Some axons within 
10 the fimbria are tracked into ventral hippocampal commissure, and axons within 
internal capsule are tracked into the corpus callosum. 

Figure 5 shows a tbree-dimensional projection and two-dimensioiuil 
validation of eight fiber bundles in the rat brain. 

With reference to Figure 5, there are shown three-dimensional 
15 projections and two-dimensional validation of 8 fiber bundles in a rat brain. The 
results of the tracking are superimposed on three-dimensional volume images using an 
oblique angle (Figure 5a) or three orthogonal angles (Figures 5b and 5c). The vertical 
Ihies between Figures 5b and 5c indicate the positions of the two-dimensional axial 
slices shown in Figures 5d-f. Color codes are; green : fimbria, dark blue : anterior 
20 commissure, light blue : medial forebrain bundle, yellow : fornix, white : stna 
terminalis^ pink : stria medullaris, red : optic tract, peach : lateral olfactory tract, A 
schematic diagram for a second embodunent of the method of the present invention is 
shown in Figure 6. In this method a fiber structure is propagated or connected in a 3D 
firom a data point of interest if adjacent pixels have anisotropy higher than a threshold 
25 value and an iimer product between the pixel of interest is higher than another 
threshold. An example of this method is shown in Figure 7. In this example, a data 
point inside an optic tract of fixed sheep brain was propagated if the adjacent data point 
had an anisotropy index of higher than 0.38 and an inner product higher than 0.98. 

The first method is suitable to delineate the detailed view of fiber 

30 structures within a white matter tract and the second method is suitable to characterize 

8 
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the global shape of each white matt^ tract. The method may also be employed with 
respect to gray matter structures. 

The present invention permits the use of standard MRI imaging 
procedures and equipment in generating the data employed to create the brain axonal 
5 fiber images of the present invention. While the invention may be employed with a 
wide range of ntiagnetic field gradient intensities, in order to maximize the signal-to- 
noise ratio, it is preferred that the magnetic field gradient be at least about 1.0 
Gauss/cm and preferably greater than 2.2 Gauss/cm 

It will be appreciated that the invention may be employed in connection 
10 with diagnostic and research uses, as well as therapeutic uses. For example, in 
connection with drug therapy employed in diseases such as Leukodystrophies 
monitoring of the fibers through the imaging procedures of the present invention will 
provide guidance with respect to ongoing therapy. 

It will be appreciated that details regarding the software employed in 
15 the present invention need not be provided in the present application as one skilled in 
the art can readily determine the software for performing the automated procedures. 
Any suitable computer, such as a naicroprocessor or personal computer, for example, 
may be employed. 

Whereas particular embodiments of the invention have been described 
20 hereinabove, for purposes of illustration, it will be evident to those skilled in the art 
that numerous variations of the details may be made without departing from the 
invention as defined in the appended claims. 



9 



wo 01/91639 PCT/USOO/14838 



CLAIMS: 

1 . A method of creating an image of brain fibers comprising 
exposing said brain fibers to a magnetic resonance imaging 

process, 

3 introducing data acquired from said magnetic resonance imaging 

process into a coicputer, 

employing a diffusion tensor to obtain a plurality of values and a 

plurality of vectors, 

employing said values and said vectors to initiate brain fiber 

10 tracking, and 

continuing or terminating said fiber tracking based upon a 
determination regarding whether the fiber is continuous or terminated based on 
randomness of fiber orientation of adjacent pixels. 

2. The method of claim 1 including 

15 employing as said values three eigen values and enq)loying as 

said vectors three eigen vectors. 

3 . The method of claim 2 including 

measurement of diffusion constant along more than six axes at 

each said pixel. 

20 4. The method of claim 3 including 

obtaining said three eigen values and three eigen vectors by 
calculating six variables in a diffusion tensor. 

5. The method of claim 4 including 

calculating said six indep^ent variables in a 3 x 3 diffusion 

25 tensor, and 

initiating said tracking at a predetennined point selected in 
three-dimensional space. 
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6. The method of claim 5 including 

effecting said tracking in orthograde and retrograde directions. 

7. The method of claim 6 including 

effecting said tracking along the eigen vector associated with the 

largest eigen value. 

8. The method of claim 1 including 

employing said conq>uter to make a determination to continue or 
terminate tracking based on the relationship 




where vxi is the unit vector representing the longest principal 
diffusion axis Xi, the index i is the three-dimensional coordinate of a pixel of interest 
(fiber coordinate) and j are the coordinate of pixels surrounding the pixel i, s is the 
number of data points referenced. In general, the fiber coordinate, i, is considered to 
have 26 surrounding points, j. A point such as I or j may be identified by three 
numbers each related to reference axis. The term "abs" refers to the number having 
an absolute value without reference to " or 

9. The method of claim I including 

employing said image information to generate a three- 
dimensional image of said brain fibers. 

10. The method of claim 8 including 

terminating said tracking when R is less than about 0.8. 

1 1 . The method of claim 1 including 

determining by said trackhig the extent and distribution of the 
projections of said fibers. 

12. The method of claim 1 including 

employing said method to create image ioformation of fibers 
disposed in the white matter portion of the brain. 

11 
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The metiiod of claim 1 including 

employing said method on a living human being. 

The method of claim 1 including 

en5)loying a continuous number field in effecting said fiber 
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